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Abstract 

The paper develops a (2 + 2)-imbedding formalism adapted to a double foliation of 
spacetime by a net of two intersecting families of lightlike hyper surf aces. The formalism is 
two-dimensionally covariant, and leads to simple, geometrically transparent and tractable 
expressions for the Einstein field equations and the Einstein-Hilbert action, and it should 
find a variety of applications. It is applied here to elucidate the structure of the character- 
istic initial- value problem of general relativity. 



1 Introduction 



The classic analysis of Arnowitt, Deser and Misner (ADM) [|IJ formulates gravitational dynamics 
in terms of the evolution of a spatial 3-geometry. The geometrical framework is the imbedding 
formalism of Gauss and Codazzi for the foliation of spacetime by spacelike hypersurfaces 0. 

Quite often, however, one encounters circumstances where a lightlike foliation is especially 
suitable. Because of the degeneracies that arise in the lightlike case the imbedding relations 
are very different and the situation not quite so familiar and under control. To bypass the 
degeneracies, one is forced to fall back to a foliation of codimension 2, by spacelike 2-surfaces. 
It is our aim in this paper to develop a simple (2 + 2)-imbedding formalism of this kind. 

Several (2 + 2)-formalisms are extant [[3], the earliest and best known being the generalized 
spin-coefficient formalism of Geroch, Held and Penrose (GHP) [||. Basically, of course, all such 
formalisms have the same content, but they take very different forms. 

The essential feature of the present approach is that it maintains manifest two-dimensional 
covariance while operating with objects having direct geometrical meaning. Two-dimensional co- 
variance permits reduction of the Einstein field equations to an especially concise and transparent 
form: the ten Ricci components are embraced in a set of just three compact, two-dimensionally 
covariant expressions. 

There is a limitation, at least in the version presented here. (It applies to most of the 
formalisms we have listed [EL) The two independent normals to an imbedded 2-surface — 
conveniently taken as a pair of lightlike vectors, since their directions are uniquely defined — are 
assumed from the beginning to be hypersurf ace-orthogonal. This precludes choosing them as 
principal null vectors of the Weyl tensor for a twisting geometry like Kerr. In this respect, the 
formalism is less flexible than GHP, and not as well tailored for the study of algebraically special 
metrics. 

(2 + 2) formalisms have a wide range of applications: to the analysis of the characteristic 
initial- value problem ||, the dynamics of strings [|J and of real and apparent horizons |7| and 
light-cone quantization || and gravitational interactions in ultra-high energy collisions 0. In a 
separate publication |L0] , we shall use the present formalism to study the nature of the singularity 
at the Cauchy horizon in a generic black hole. 

We conclude this Introduction by briefly outlining the contents of the paper. The basic 
metrical notions (adapted co-ordinates, basis vectors and form of the metric) are defined in 
Sec. 2. In Sec. 3, we introduce in two-dimensionally covariant form the geometrical information 
encoded in first derivatives of the metric: the extrinsic curvatures and "twist," as well as the 
invariant operators which perform differentiation along the two lightlike normals. This comprises 
the basic formal machinery needed in Sec. 4, which presents the central result of the paper, the 
tetrad components of the Ricci tensor as three concise equations (p7|) (p9|) . (To make direct 
access to these results easy, their derivation is deferred to the second half of the paper (Sees. 9- 
12), which also provides (Sec. 13) the tetrad components of the full Riemann tensor.) 

The contracted Bianchi identities (Sec. 5) are applied in Sec. 7 to analyze the structure of the 
characteristic initial-value problem. In Sec. 8 we sketch the Lagrangian formulation of covariant 
double-null dynamics. 

The Ricci and Riemann components result from the commutation relations for four-dimen- 
sional covariant differentiation. Their most efficient derivation calls for a formalism that is both 
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four- and two-dimensionally covariant. Unfortunately, these two requirements do not mesh easily. 
Four-dimensional covariance tends to clutter the formulae by treating subsidiary two-dimensional 
quantities like shift vectors and the two-dimensional connection as 4-scalars, on a par with the 
primary geometrical properties, extrinsic curvature and twist. Those properties, for their part, 
are correlated, not with four-dimensional covariant derivatives, but with Lie derivatives, which 
are non-metric and have no direct link to curvature. To patch up these differences, and thus 
streamline the derivations, seems to need a certain degree of artifice. In Sec. 10 we address this 
(purely technical) problem by temporarily working with a "rationalized" covariant derivative 
which exhibits both four-dimensional and restricted ("rigid") two-dimensional covariance. 
Some brief remarks (Sec. 14) conclude the paper. 



2 (2 + 2)-split of the metric 

We shall suppose that we are given a foliation of spacetime by lightlike hypersurfaces S° with 
normal generators ££>, and a second, independent foliation by lightlike hypersurfaces X 1 with 
generators tf£ nowhere parallel to £^'. The intersections of {S } and {S 1 } define a foliation of 
codimension 2 by spacelike 2-surfaces S. (The topology of S is unspecified. All our considerations 
are local.) S has exactly two lightlike normals at each of its points, co-directed with £^ and 
£W. 

In terms of local charts, the foliation is described by the imbedding relations 

x a = x a (u A ,9 a ). (1) 

Here, x a are four- dimensional spacetime co-ordinates (assumed admissible in the sense of Lich- 
nerowicz |TT[); u° and u 1 are a pair of scalar fields constant over each of the hypersurfaces S° 
and S 1 respectively; and 8 2 , 6 3 are intrinsic co-ordinates of the 2-spaces S, each characterized 
by a fixed pair of values (u°, u l ). 

Notation: Our conventions are: Greek indices a, j3, . . . run from to 3; upper-case Latin 
indices A,B,... take values (0, 1); and lower-case Latin indices a,b, . . . take values (2, 3). We 
adopt MTW curvature conventions with signature ( — h ++) for the spacetime metric g a p. 
When there is no risk of confusion we shall often omit the Greek indices on 4-vectors like £^ 
and e? a \ : they are easily identifiable as 4-vectors by their parenthesized labels. Four- dimensional 
covariant differentiation is indicated either by V a or a vertical stroke: VfsA a = A a \p. Four- 
dimensional scalar products are often indicated by a dot: thus, £^a) '^(b) = 9a.p£?A) ^(b)- Further 
conventions will be introduced as the need arises. 

Without essential loss of generality we may assume the functions x a (u ,8 a ) to be smooth 
(at least thrice differentiable). (We are always free to make the co-ordinate choice x A = u A , 
x a = 8 a , but at the cost of losing manifest four- dimensional and two-dimensional covariance.) 

The lightlike character of the hypersurfaces T, A is encoded in 

V/ ■ Vu B = g a/3 (d a u A )(dpu B ) = e-\ AB (2) 
for some scalar field X(x a ), where 

rj AB = anti-diag (-1, -1) = r] AB ; (3) 
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n AB and its inverse 7]ab are employed to raise and lower upper-case Latin indices, e.g., £( ) = 
-£('). 

The generators l^ A ' of S A are parallel to the gradients of u A (x a ). It is symmetrical and 
convenient to define = e x Vu A , i.e., 

4 A ) = e x d a u A . (4) 

Then 



£ (A) • = eH B A . (5) 
The pair of vectors er a ), defined from ([!]) by 

ef a) = dx a /d6 a , (6) 

are holonomic basis vectors tangent to S. The intrinsic metric g a b d6 a d9 b of 5* is determined by 
their scalar products: 

9ab = e (a) ■ e (6) . (7) 

Lower-case Latin indices are lowered and raised with g a ^ and its inverse g ab ; thus = g ah C(b) 
are the dual basis vectors tangent to S, with • en,) = 5%. Since £^ A ^ is normal to every vector 
in T, A , we have 

i {A) ■ e (a) = 0. (8) 




Figure 1: The 2 + 2 splitting of the four dimensional space time into a foliation of 
intersecting null surfaces u° = const, and u 1 = const.. 

In general, 9 a cannot be chosen so as to remain constant along both sets of generators £^ A \ 
They are convected (Lie-transported) along the pair of vector fields dx a /du A (in general, non- 
lightlike). 
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From (JD and (J5]) one finds that dx a /du A — l? A s is orthogonal to £(b), i.e., tangent to S. This 
validates the decomposition 

dr a 

thus defining a pair of "shift vectors" s A tangent to S (see Fig. 1). 

An arbitrary displacement dx a in spacetime is, according to ([|) and (|J), decomposable as 

dx a = t a (A) du A + ef a) (d6 a + s a A du A ). (10) 

From (||), (0) and (|8|) we read off the completeness relation 

fin 



Combining (fL0|) and ([TT|) shows that the spacetime metric is decomposable as 

g a/3 dx a dx p = e x r)AB du A du B + g ab (d6 a + s a A du A )(d6 b + s b B du B ). (12) 

3 Two-dimensionally covariant objects embodying first 
derivatives of the metric: extrinsic curvatures KAab-, 
twist u a and normal Lie derivatives Da 

Absolute derivatives of four-dimensional tensor fields with respect to u A and 6 a are projections 
of the four- dimensional covariant derivative V«, and denoted by 

5 dx a „ 5 „ 

V «> TTTn = e (a) ■ V. (13) 



5u A du A °' 56 a 
From d^) and the symmetry of the mixed partial derivatives and the affine connection 

5u A ~ 56 a \du A ) ' 56 b ~ 56 a ' 

The object 



(14) 



T c ah = e^-5e {a) /56 b (15) 

is, as the notation suggests, the Christoffel symbol associated with g a b, as is easily verified by 
forming d c g a b, recalling (|7|) and applying Leibnitz's rule. 

Associated with its two normals £(a), S has two extrinsic curvatures K^ab, defined by 

K A ab = e (a ) ■ 5£ {A) /56 b = t(A) a \p ef a) e p {b) . (16) 

(Since we are free to rescale the null vectors it a), a certain scale-arbitrariness is inherent in this 
definition.) Because of @, we can rewrite 

K Aab = -1 {A) ■ 5e {a) /56 b , (17) 

which exhibits the symmetry in a, b. 
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A further basic geometrical property of the double foliation is given by the Lie bracket of 
i(B) and £(A), i-e., the 4- vector 



[i (B) ,£ {A) } a = 2(£ m .V)£f A)] . (18) 

Noting (|S|) and the fact that the Lie bracket of the vectors dx a /du B and dx a /du A vanishes 
identically, and recalling (0), we find 

[£(B)J(A)] = e A BUJ a e {a)} (19) 

where the 2-vector uo a is given by 

u a = e AB (d B s A -s B s% b ), (20) 

the semicolon indicates two-dimensional covariant differentiation associated with metric g a b, and 
€ab is the two-dimensional permutation symbol, with eoi = +1. (Note that raising indices with 
rj AB to form e AB yields e 10 = +1.) 

The geometrical significance of the "twist" u a can be read off from (19): the curves tangent 
to the generators £(p) , £(i) mesh together to form 2-surfaces (orthogonal to the surfaces S) if and 
only if u a = 0. In this case, it would be consistent to allow the co-ordinates 9 a to be dragged 
along both sets of generators, and thus to gauge both shift vectors to zero. 

We denote by D A the two-dimensionally invariant operator associated with differentiation 
along the normal direction £i A y Acting on any two-dimensional geometrical object X a b " , D A is 
formally defined by 

d a x\ : . = (o A - c s , A )x a b :_. (21) 

Here, d A is the partial derivative with respect to u A and C s d the Lie derivative with respect to 
the 2-vector s d A . 

As examples of fl2"T|), we have for a 2-scalar / (this includes any object bearing upper-case, 
but no lower-case, Latin indices): 

D A f = (d A - s a A d a )f = £ a {A) d a f (22) 

(in which the second equality follows at once from (P)); and for the 2-metric g a b'. 

D A g a b = d A g a b — 2s A ( a -b) = 2K Aab , (23) 

in which the second equality is derivable from (|7|), (|i"4|) and (For the detailed derivation, see 
(ff5|) below, or Appendix B.) 

The geometrical meaning of D A is quite generally the following (see Appendix B): 
D A X a b is the projection onto S of the Lie derivative of the equivalent tangential 4-tensor 

ya... — ya... a (6) 
A 13... = A b... e (a) e l3 ■ ■ ■ 

with respect to the 4-vector if 



'{Ay 



D A X a bZ = e^e? b) ■ Cq Xfz ■ (24) 
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The objects K Aa b, u a and Da comprise all the geometrical structure that is needed for a suc- 
cinct two-dimensionally covariant expression of the Riemann and Ricci curvatures of spacetime. 
According to (|IED , (|T9|) and fl24|), all are simple projections onto 5* of four- dimensional geo- 
metrical objects. Consequently, they transform very simply under two-dimensional co-ordinate 
transformations. Under the arbitrary reparametrization 

6 a -> 6 a ' = f a (6\u A ) (25) 

(which leaves u A and hence the surfaces E" 4 and S unchanged), u a and K^ab transform cogredi- 
ently with 

e (a) - e[ a) = e (b) d6 b /d6 a ' (26) 

(see (|6|)), and Da is invariant. By contrast, dx a /du A and hence the shift vectors s a A (see ([|)) 
undergo a more complicated gauge-like transformation, arising from the ^-dependence in fl2~5"|). 

4 Ricci tensor 

The geometrical and notational groundwork laid in the previous sections allows us now to simply 
display the components of the Ricci tensor, deferring derivations to Sees. 9-12. Our notation for 
the tetrad components is typified by 



^Rab — Rape'(a) e (b)i ^aA ~ Rafi^a)^ 



(A)- 



The results are 



MR* = l -^Rg ab -e- x {DA + K A )K A 

+2e- x K A ^ a K b A )d - l -e- 2X u a u h - \, ab - l -\ a \ b (27) 

RAB = —D^aKb) — KAabKs ab + K(aD B )\ 



R Aa = K2 a;b -d a K A --d a D A X + -K A d a X 

1 A 

+ -e AB e- x 



D + K*) D E \ - e- A uo a u a + (e A ) ;a a (28) 
1„ „ . 1 



[D B + K B ) uo a - uo a D B \\ , (29) 
where ^R is the curvature scalar associated with the 2- metric g ab , and Ka = K / 



L A a- 



5 Bianchi identities. Bondi's lemma 

The Ricci components are linked by four differential identities, the contracted Bianchi identities 

VpR^ = h a R, (30) 
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where the four-dimensional curvature scalar R = R a a is given by 

R = e~ x R A A + R a a , (31) 

according to fllTD- 

As we show in Sec. 12, projecting (|30|) onto e( a ) leads to 

(D A + K A ) R A a = l -d a R A A + l -e x d a WR\ - (e x ^R b a ) & . (32) 
Projection of fl3"0|) onto similarly yields 

Db (R\-\^R D D )-\e X D A ^R\ 

= e x ^R ab Kf - R\K B - (e x R\) ^ + e AB u a R B a . (33) 



Equations (p2|) and (|33D express the four Bianchi identities in terms of the tetrad components 
of the Ricci tensor. 

We now look at the general structure of these equations. 

For A = in (^), R® does not contribute to the first (parenthesized) term, since 

— Roi = Rq = R\ = 2^a- (34) 
This equation therefore takes the form 

D.Roo + l -e x D ^Rl = -K,R m + C(^R ab , i? o, R 0a , d a ), (35) 

in which the schematic notation C implies that the expression is linear homogeneous in the 
indicated Ricci components and their two-dimensional spatial derivatives d a . 
The other (A = 1) component of ( |3~3"D has the analogous structure 

D Q Ru + ~e A A (4) K = -K x Roi + £( ( %„ 6 , R n , Ria, da)- (36) 
The form of the remaining two Bianchi identities (|32|) is 

DoRla + DlRoa = £(} Rab, Roi, RAa, 9 a ) . (37) 



It is noteworthy that the appearance of i?oi i n ( p5|) and (36) is purely algebraic: its vanishing 
would be a direct consequence of the vanishing of just six of the other components. Bondi et al 
13|] and Sachs therefore refer to the _R i field equation as the "trivial equation." 

The structure of (|35"l)-(j37l) provides insight into how the field equations propagate initial data 



given on a lightlike hypersurface. Let us (arbitrarily) single out u° as "time," and suppose that 
the six "evolutionary" vacuum equations 

^R ab = 0, R 00 = R 0a = (38) 
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are satisfied everywhere in the neighbourhood of a hypersurface u° = const. (Bondi and Sachs 
refer to ^R a b as the "main equations" and to Rqo, Ro a as "hypersurface equations." Rqo is, in 
fact, the Raychaudhuri focusing equation |12|], governing the expansion of the lightlike normal 
£( ) = — ft- 1 ' to the transverse hypersurface u 1 = const., and Ro a similarly governs its shear.) 

Then ( 55 ) shows that the trivial equation R m = is satisfied automatically. From (|5E|) and 
fl37|) it can be further inferred that if i?n and R\ a vanish on one hypersurface u° = const., then 
they will vanish everywhere as a consequence of the six evolutionary equations (§^). This is the 
content of the Bondi-Sachs lemma |TJ|, [J, which identifies the three conditions R\\ = R\ a = 
as constraints — on the expansion and shear of the generators £m of an initial hypersurface 
u° = const. — which are respected by the evolution. 



6 Co-ordinate conditions and gauge-fixing 

The characteristic initial-value problem |j) involves specifying initial data on a given pair of 
lightlike hypersurfaces E°, E 1 intersecting in a 2-surface Sq. 

It is natural to choose our parameters u A so that u° = on E° and u l = on E 1 . The require- 
ment (0) that u A be globally lightlike already imposes two co-ordinate conditions on (u A ,9 a ), 
considered as co-ordinates of spacetime. Two further global conditions may be imposed. We may, 
for instance, demand that 8 a be convected (Lie-propagated) along the lightlike curves tangent 
to £f Q \ from values assigned arbitrarily on S°. According to @, this means the corresponding 
shift vector is zero everywhere: 

S( ) = -ifa dJ a = 0. (39) 

In this case, (^0) shows that 

u a = d s1 (40) 

is just the "time" -derivative of the single remaining shift vector. 

These global co-ordinate conditions can still be supplemented by appropriate initial condi- 
tions. We are still free to require that 6 a be convected along generators of E° from assigned 
values on So] then 

4 = (u° = 0) (41) 



in addition to (p9|). 

In addition to (or independently of) (139|) and (fUD, we are free to choose u 1 along S° and u° 
along E 1 to be affine parameters of their generators. On E°, for instance, this means, by virtue 
of (|) and®, 

gen. 



so that A vanishes over S°. There is a similar argument for E 1 . Thus, we can arrange 

A = (E° and E 1 ). (42) 
Alternatively, in place of (pffi), the co-ordinate condition 

D x \ = -K x on E° (43) 
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could be imposed to normalize u l . (A corresponding condition on S 1 would normalize u°.) 
The Raychaudhuri equation (p8|) for Ru on S° would then become linear in the expansion rate 
K\ = d\ In (72, and that facilitates its integration (cf Hayward ||, Brady and Chambers J/J). 



7 Characteristic initial-value problem 

We are now ready to address the question of what initial data are needed to prescribe a unique 
vacuum solution of the Einstein equations in a neighbourhood of two lightlike hypersurfaces S° 
and S 1 intersecting in a 2-surface S 0. 

We arbitrarily designate u° as "time," and shall refer to S° (u° = 0) as the "initial" hyper- 
surface and to S 1 {u l = 0) as the "boundary" 

We impose the co-ordinate conditions (H), §3) and g|) to tie down 6 a and u A . While fl3|) 
and (|4l|) control the way 6 a are carried off S , onto S° and into spacetime, the choice of 6 a on So 
itself is unrestricted. Thus, our procedure retains covariance under the group of two-dimensional 
transformations 9 a — > 9 a = f a (6 h ). 

In the 4- metric g a/ 3, given by (|I2D, the following six functions of four variables are then left 
undetermined: 

9ab, A, s a v (44) 

(In place of sf, it is completely equivalent to specify u a = dos®, since the "initial" value of si is 
pegged by (0).) 

We shall formally verify that a vacuum 4-metric is uniquely determined by the following 
initial data: 

(a) . On S*o, seven functions of two variables 6 a : 

g ab , u a , K A = d A hyg$ (S ); (45) 

(b) . on S° and S 1 , two independent functions of three variables which specify the intrinsic 

conformal 2- metric: 

g-iga, (S° and S 1 ). (46) 



Instead of (f4q ), it is equivalent to give the shear rates of the respective generators, 

^onS°, ^onS 1 , (47) 
defined as the trace-free extrinsic curvatures: 

&Aab = K Aab - ~ g a bK A = )^g^d A (g^ g ab ). (48) 

These two functions correspond to the physical degrees of freedom ("radiation modes") of the 
gravitational field || |i"4| . 

To build a vacuum solution from the initial data ([45]), (f47D, we begin by noting that (|39| ) 
implies that D = do, K 0a b = K 0ab everywhere. Hence the general expression (|28D for R 00 
reduces here to 

- R 00 = (d + \Ko - \o) Ko + aoabcrf- (49) 
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On E 1 , we have A = A o = by fl42|). Thus, ([19|) becomes an ordinary differential equation for 

K = K = d \iag^ 



as a function of u°. This can be integrated along the generators, using the given data for ao b a on 
E 1 , and the initial value of K on Sq, to obtain gz , hence the full 2-metric g a b (hence also K§ a b) 
over E 1 . 

Expression (^9|) for i? 0a reduces similarly to 



Roa = ~\e' x (d + K - A, ) u a -\{d Q - K )\ a + K^ b - d a K (50) 

in a spacetime neighbourhood of E° and E 1 . On E 1 , since K 0ab is now known, and A = A i0 = 
A ; o = 0, ( |50| ) is a linear ordinary differential equation for uj a which may be integrated along 
generators, with initial condition (45), to find u a (hence si). 



Thus, our knowledge of the six metric functions ( pH|) has been extended to all of E 1 with the 
aid of the evolutionary equations R 00 = R 0a = 0. 

A similar procedure, applied to the constraint equations Rn = R\ a = 0, determines the 
functions (|44]) (hence also K\ a b) over the initial hypersurface E°. (Here we exploit (f41~D — implying 
D\ = d\ — which holds on E° only. This limitation is of little practical consequence, since the 
Bianchi identities (Sec. 7) relieve us of the need to recheck the constraints off E°.) 

Thus, the data flUf) , together with their tangential derivatives d\, d a — which we denote 
shematically by 

T> = {g a b, A, u a , si, d u d a } (51) 



-are now known all over the initial hypersurface E u = 0. (Note that T> includes K 



lab- 



We now proceed recursively. Suppose that T> is known over some hypersurface E : u° = const. 
We show that the six evolutionary equations ^R a b = 0, -Roo — Roa = 0, together with the known 
boundary values of g a b, K 0ab , uj a and on E 1 , determine all first-order time-derivatives d of V, 
and hence the complete evolution of T>. 

Expression (27) for the evolutionary equations ^R a b = can be written more explicitly, with 



the aid of the identity 

D A K* - 2K A ( a K b A )d = -2D 1 K 0ab + AK 0( d a K lb)d + u (a]bh (52) 
which is rooted in the symmetry 

9[BK A ] ab = 0, K Aab = K Aab + SA(a;fe) (53) 

(see (p3|) and Appendix B). 

The equations ^R a b — are thus seen to reduce to a system of three linear ordinary differen- 
tial equations for K 0ab as functions of u 1 on E, whose coefficients are concomitants of the known 
data T> on E. Together with the boundary conditions on K 0ab at u 1 = (i.e., the intersection of 
E with E 1 ), they determine a unique solution for K 0ab on E. 

We next turn to (|49|) and ( |50D to read off the values of <9 A and d u! a on E. Since the remaining 
time-derivatives are known trivially from 

d s1 = u a , ^d g ab = K 0ab = K 0ab , 
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we are now in possession of the first time-derivatives of all the data DonS. 

This completes our formal demonstration that the initial conditions fl45|) and (fR|), or (|45| ) 
and (pETD, determine (at least locally) a unique vacuum spacetime. 



8 Lagrangian 

According to ( |12|) and (|31~1) , the Einstein-Hilbert Lagrangian density £ = (— A g)^R^ decomposes 

£ = gh\e- x Ri+ ^R a a ), (54) 



as 



in which <p refers to the determinant of g a b- Substitution from ( f27|) and (p8|) yields the explicit 
form 



g-hc = e x (2) i? - D A {2K A + D A \) - K A K A - K a A b K A b 



e~ x cu a tu a -e x 2A ;a a + -A, a A-° 



(55) 



Second derivatives of the metric in (|55|) can be isolated in the form of a pure divergence by 
calling on the identities 



g*D A x j 



-*g)le- x X A < 



'{A) 



- g l *X A K A , 



A a , a + A a \ a = V Q (A a ef a) ), 
which follow from ( |102| ) below, and hold for any scalars X A and 2- vector A a . We thus obtain 



(56) 
(57) 



£ = -d a 



+9 l 



a 
(a) 



(- *g) 5 e x (2K A + D A \) t {A) + 2 (- *g) 1 A< a e 
e x + K A K A - KfK A + V A c^ a + K A D A \ + ^e A A, a A' a 



(5* 



The divergence term integrates as usual to a surface term in the action S = J £ d A x, and has no 
influence on the classical equations of motion. 
Variation of S with respect to 

-e A = g(oxi) = £( ) • £{i) 

reproduces the expression obtained from (P7|) for G 01 = ^e x R^. Similarly, variation with respect 



to s A yields the expression (P9|) for G Aa = R Aa , if we take account of the implicit dependence of 
K Aab , D A and uj a on s a A through 



K Aa b = ~d A g a b — SA(o;6)> 

2|) and Finally, variation with respect to g a j, yields ( 4 )G ab , if we note the identity 



(59) 



ab<P' c~ (P;ab- 



(60) 
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Thus, variation of the action fl5"8|) yields eight of the ten Einstein equations. The remaining 
two equations — the Raychaudhuri equations for i? o and Ru — cannot be retrieved directly from 
(pBD , because the a priori conditions rj 00 = rj 11 = (expressing the lightlike character of u° and 
u l ) which is built into (|58p, precludes us from varying with respect to these "variables." The two 
Raychaudhuri equations can, however, be effectively recovered from the other eight equations 
via the Bianchi idendities. 

The Hamiltonian formulation of the dynamics has been discussed in detail by Torre ||. We 
hope to pursue this topic elsewhere. 

9 Gauss- Weingarten (first order) relations 

In this second half of the paper, we return to the beginning and to the task of laying a more 
complete geometrical foundation for the Ricci and Bianchi formulas which we quoted without 
derivation in (p7|) (p9|) and (0), (|33|). We begin with the first-order imbedding relations for the 
2-surface S subspace of spacetime. 

(|T5|) and ( |T7j ) allow us to decompose the 4- vector Se? a JS9 b in terms of the basis (£(A),e(a))- 



Recalling (g), we find 



S9 b 

Similarly, in view of ([16]), we may decompose 



--e~ X KaAA)+K b e {c) . (61) 



^ = L AB a £W + K Aab e (b) (62) 

where the first coefficient is given by 

Lab* = e~H {B) ■ 8£ {A) /S9 a . (63) 
This coefficient can be reduced to a much simpler form. Its symmetric part is 

L(AB)a = ^ e~ X d a {?(A) • £&)) = ^n AB d a \. (64) 
To obtain the skew part, we note first that 

e W Z[(B)t(A)}la\beta] = - (65) 



since £( A ) is proportional to a lightlike gradient (see (^)). With the aid of (|65|) and (^9|) we now 
easily derive 

L[AB]a = e- x ^ {B) £ {A )]i3\aef a) = -e~ x [£ {B) ,£ { A)} a e {a)a 

= ^e A B^ a e~ x , (66) 
Combining (|64]) and (§3), we arrive at the simple expression 

2L ABa = rj AB da^ + e~ x t AB u a (67) 
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for the first coefficient in (|6^). 

The Gauss- Weingarten equations ( |6lD and (|62|) govern the variation of the 4- vectors £( A ) and 
e( a ) along directions tangent to S. We now turn to their variation along the two normals. 

We have from (Q), 

V p i [A)a = 2£ {A)[a d f3] X + V a £ {A)/3 . (68) 
Multiplying by £^ B ^ and symmetrizing in A, B gives 

(t(B) ■ V)^ ( A) + (£ { A) ■ V)£ (B) = 2£ ({A) D B) X - r] AB e x VX (69) 

where D A \ is defined as in (|22|) . It follows that 

VA = e~ x £ {A) D A X + e (a) d a X. (70) 

On the other hand, the difference of the two terms on the left of ( |69"D is given by ( |19|) as e AB u a e^ . 
Adding finally yields 

(£ {B )-V)£ { A) = N ABC £ {C) -e x L BA J a \ (71) 

where ^ 

N ABC = D {A Xr] B)c - -r] AB D c X, (72) 

and L was defined in fl5"Tp. 

Proceeding finally to the transverse variation of e( a ), we have from (|14]) and (|9]), 

Substituting from (^) and (|62|) , it is straightforward to reduce this to 

^ = ( L ABa - s b A e- x K Bab )£^ + K Aa b e (6) (73) 

where 

K Aab = K Aab + s^.o. (74) 

Applying Leibnitz's rule to 
we read off from 0731) the result 



9 A9ab = JffX^) ' e (6))> 



-^d A g ab = K Aab = K A (a b ), (75) 

which gives direct geometrical meaning to the extrinsic curvature in terms of transverse variation 
of the 2- metric. 

The normal absolute derivatives of e( a ) are given by (recalling @ and (|T3|) ) 

(Aa) " V)e (a) - - s A - w . 
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With the help of (|73| ) and ([JTj) this reduces to 

(i(A) ■ V)e (a) = L AB + (K Aa b - s A c T h ac )e [b) (76) 

Correspondingly, the two normal derivatives of g ab are 

(£(A) ■ ^)9ab = 2K Aab - s c A d c g ab . (77) 

The two-dimensionally noncovariant terms which appear in ( [76] ) and ( |77| ) are not a mistake. 
They arise because the normal gradient £( A \ ■ V, applied to objects carrying lower-case Latin 
indices — let us say g ab — does not preserve manifest two-dimensional covariance, since it contains 



(see (|22|)) a piece —s c A d c g ab involving ordinary (rather than two-dimensional covariant) derivatives 
with respect to 9 C . Although not incorrect, this is a formal impediment: it threatens to clutter 
our formulae with terms in the shift vectors s a A which are, to boot, noncovariant. In the following 
section, we explain how this can be remedied by introducing a "rationalized" gradient operator 
V. 



10 Rationalized operators V, Da, V a 

The rationalized operator V Q avoids the two-dimensionally noncovariant terms which appear 
when V a is applied to objects bearing lower-case Latin indices, as in (plj), ( [75] ) and fl77|). 

Applied to scalar fields or to 4-tensors not bearing lower-case Latin indices, V is identical 
with V. If the object does carry such indices, there are supplementary terms involving the 
two-dimensional connection T£ c . 

Specifically, we define 

V a = V a +pi a) (Va-e (a) -V) (78) 

in which eu) ■ V = 5/56 a is the absolute derivative introduced in ([TBI), and the operator V a will 
be specified in a moment. We have introduced the pair of 4-vectors = V# a , i.e., 

p (a) = dQa/Q x *. (79) 

Their projections onto the basis vectors are, according to (||) and 

p<°> • e (6) = 81 p M . £ {A) = - s % (80) 
from which follows the identity 

€-P^eU = ^i A) d^/du A . (81) 
Hence ([78] ) can be recast in terms of the absolute derivative S/Su A : 

V = e- x £ {A) d/du A + p {a) V a . (82) 
We next introduce the differential operator 

D A = £ {A) ■ V = e (A) ■ V - 4(V a - 5/56 a ). (83) 
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An alternative form 

D A = 6/5u A - s a A V a (84) 



follows from fl82). 

Since e( ) • V = V a , we can reconstruct V from ( |83"D in yet another form: 

V = e- x e {A) D A + e (a) V a . (85) 

We now specify the operator V a . It is defined so as to act as a two-dimensional covariant 
derivative on all lower-case Latin indices (including parenthesized ones), and at the same time 
as an absolute derivative 5/59 a on Greek indices. Upper-case Latin indices are treated as inert. 

As an example, 

V b e {a) = 5e (a) /56 b -Tl b e {c) . (86) 

It is evident that, quite generally, the "correction" V a — 5/59 a in ( |83| ) and (|7|) is linear and 
homogeneous in the two-dimensional connection T^ b . 

Examples of how V a and D A act on scalars and 2-tensors are 



Vf) 



VJ = d a f, D A f = (d A - s a A d a )f = D A f, 

v ( ,v , : :: . = Xt, a , D D A X h Z _ = {d A - s\V a )X h Z _. 

For the 2-metric g ab , we have from d8~4|). 

D A g ab = d A g ab , (88) 
since V ' c g a b = 9ab;c = 0. Thus, (|75|) can be expressed as 

]^D A g ab = K Aab , (89) 

which should be contrasted with (|77|) . 

Similarly, with the aid of (|86| ) and flSB]), the noncovariant expressions ( |6TD and (^) become 

V 6 e (a) = V a e (b) = -e- x K Aab £ (A \ (90) 

D A e {a) = L ABa i {B) + K Aab e {b) . (91) 

Quite generally, D A , V a and V preserve both four-dimensional covariance and covariance un- 
der "rigid" two-dimensional co-ordinate transformations 9 a —>■ 9 a = f a (9 b ), with no dependence 
on u A . (w-dependence of f a would induce "gauge" transformations of the shift vectors s A , see 
the remarks following fl26D.) 

With the aid of (|85D , the last two results can be put together to form the rationalized covariant 
derivative of e^ a y. 

e x V p e {a)a = {L BAa t^ + K Bab e^f ] - K Aah 6 A) ef. (92) 

Equations ( |62j) and (|7T|) similarly combine to produce 

V^ (A) a = (e- X N ABC i£) - L BAa e^)if + (L ABb £^ + K Aab e^)ef . (93) 

(No distinction here between V and V, since ti A ) carries no lower-case Latin indices.) 

To sum up: equations (|92|) and (^) encapsulate the full set of first-order ( Gauss- Weingar- 
ten) equations, which control tangential and normal variations of the basis vectors e( Q ), £r A y The 
coefficients in these equations are given by (|16"1) , (|67|), (|20|), ([72] ) and ([741). Their geometrical 
meaning emerges from ([75]), ( |i~9]) and remarks following those equations. 
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11 Rationalized Ricci commutation rules 

The usual commutation relations need to be modified for V a . To derive the modified form, 
consider the action of V on any field object X a bearing just one lower-case Latin and an arbitrary 
set of other indices. 
From (|7§) and (§§), 

v 7 v^x a = (5y, - p^K a )(8 b yp - pf ] v b mc )x b . 

Skew-symmetrizing with respect to (3 and 7, and noting from (|79|) that V^p^ = leads to 

V h V p] X a = V h V^X a +p^ ) {p^R b amn - e- x dtf d A T b ma }X b , (94) 

in which <9 7 r^ a has been expanded using 

d^ = e- x ^d A +p^d n , (95) 

which is a special case of (|82|) . 

The right-hand side of (Q) can be further reduced: d A Y b ma is a 2-tensor, given by 

d A^ b ma = 2K \a;m) ~ K Ama ' & (96) 

according to (|75|) ; and in two dimensions we have 

R b amn = R& b \ m 9n]a- (97) 



If e? a s is substituted for X a , (0), fl92|) and C|93f) can be used to express the projection onto e( Q ) 
of the four- dimensional Riemann tensor in terms of the first-order Gauss- Weingarten variables 
K, L, N and their derivatives. If our interest is primarily in the Ricci tensor, the contracted 
form (7 = a) of fl9~l| ) suffices: 

(V«V„ - V^Vje^ = ef a) R aP - i ®Rp {a)p + e- x (d a K A )£ ( p A) (98) 

where 

K A = K a Aa = d A \ng^ (99) 

and g = det g a b- 



12 Contracted Gauss-Codazzi (second order) relations. 
Ricci tensor 

The Gauss-Codazzi relations are the integrability conditions of the system of first order (Gauss- 
Weingarten) differential equations fl92|), (p3|). As just noted, they express projections of the 
four-dimensional Riemann tensor in terms of K, L, N and their first derivatives. 
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Contraction of these equations gives frame components of the Ricci tensor. The most con- 
cise way of deriving these components in practice is through recourse to a generalized form of 



Raychaudhuri's equation [12]. 



Let A a be an arbitrary 4-vector (which may bear arbitrary label indices) and B a a second 
vector free of lower-case Latin indices, so that V ' pB a = VpB a and the standard commutation 
rules apply. Then it is easy to check the identity 

R aP A a B p = V p {A a V a B?) - A a V a (VpBP) - {V p A a ){V a B^). (100) 

If, on the other hand, B is replaced by en,), then we call upon the commutation law (^3) for 
V, with the result 

R a pA a e[ b) = Vp{A a V a e P {b) ) - A«V a {V^ (b) ) - (V^ a )(V Q ef fe) ) 

+~ ®R(A ■ p {b) ) - e-\d b K B )(A ■ (101) 

With the choices A = £/ A ) and et a ), B = tt B ) we can recover all frame components of the 
Ricci tensor from these equations in tandem with (|92|) and (|93]) . 

Some details of these calculations are recorded in Appendix A. The final results have already 
been listed in (0)-©. 

We next turn to the contracted Bianchi identities. The projection of ( j30|) onto et a ) yields 

W p {R^ {a) )-RiV p e a {a) = \d a R. 
The second term is evaluated with the aid of (0). In the first term, we expand 

and note the (often used) results 

V p ef b) = d b X, V p f {A) =K A + D A X, (102) 
which follow from ( p2|) and (|93|) . The result is (p2|), and ( p3|) is obtained similarly. 

13 Riemann tensor 

We list here the tetrad components of the Riemann tensor, obtainable from the uncontracted 
Ricci commutation rules (see, e.g., (p4|)). The notation for the tetrad components is as in Sec. 4. 

^R ab cd = MR5^ b d] -2e- x K/ [c K A % 

Rabcd = ^e AB e CD (2e x D E D E X-3uj a u a + e 2X X' a X, a ) 

RAabc = 2K Aa [ b]C \ ~ K Aa [ b X jC ] ~ e~ X € AB K E b UJ c ] 
RaABC = \ e BC {D A UJa + K Aab u b - e x e AE {D E d a X - K E X' b ) - U a D A X] 

R\ B b = -d^kS + k^ + d^xkS-Iv^DeXK^ 

1 1 1 

--V AB (e' X uJaUJ b + e A A, Q A, fe + 2e A A ;ab ) - -e AE ^ e ab r. 
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We have here defined 

e ~\ = g-$ e ab d a {e~ x uj h ) 

14 Concluding remarks 

The (2 + 2) double- null imbedding formalism developed in this paper leads to simple and geomet- 
rically transparent expressions for the Einstein field equations (|2"7l)-(|2"9D and the Einstein-Hilbert 
action fl96|). It should find ready application in a variety of £1X6 clS, 8lS indicated in the Introduction. 

(2 + 2) formalisms are certainly not new || |J, but they have languished on the relativist's 
back-burner. We hope that this exposition will play a role in promoting these verstile methods 
from the realm of esoterica into an everyday working tool. 
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A Computing Ricci components: some intermediate de- 
tails 



For the convenience of enterprising readers who wish to derive the Ricci components (^4| 
for themselves, we record here some intermediate steps of the computations. 
Computation of ^R a b from ( |1U1| ) requires evaluation of 

V^(e? a) V Q e? M ) = -e-\D A + K A )K* 



(V^f a) )(V Q ef 6) ) = ^{\ a \ b + e- 2X u a u b )-2e- x K A ^ a K A b)d) 
which can be verified from (EO) and (pTp. 



Computation of Rab from ( |100|) requires 

I 

2 



«(Wf|) = (D A X)(D B X)--r, AB (D E X)(D E X) 



+K Aab K a B b - ^e x r] AB (X ta X' a + e^u^) 



which follows from (0), (0) and 

Finally, computation of R Aa requires 



(V/3 ef a) )f {A)]a = ~ \ a D A X + (K Aab + \AK Aab 



2 ' V 2 

+ l -t AB e-\uj a D B \ + AK*cu b ) 

in which 

AK Aab = K Aab — K Aab = s Ab - a . 
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B The operator Da> commutation rules and other prop- 
erties 

In Sec. 3 we gave two definitions — (pip and (|24l) — for the operator D A . It is straightforward to 
show their equivalence. We have 

[e(a),e (6) ] = 0, [dx/du A , e (6) ] = 0, 



since the Lie bracket of two holonomic vectors vanishes (cf (|14|) ). In combination with (|9]), this 
yields 

[£(A), e( 6 )] = -[s A e {a ), e {b) } = (d b s a A )e (a) . (Bl) 

Hence, for any 2- vector X b , 

C kA) {X b e {b) ) = {(d A - C s a A )X b } e (6) (B2) 

which proves the equivalence of fl2"4|) and (|2l|) when applied to X b . 

This argument is easily extended. For instance, for the 2-metric g ab , the definition fl2~j|) gives 

D A g ab = ef a) e^Ct {A) ga/} = 1e a {a) e (6) ^(A)(a|/3) 
= 2K Aab 



by (0), which agrees with the form (^) obtained from the definition ([H]). (Strictly speaking, 
( j24|) requires that the projector A aj3 = e^ a should replace g a p in the first equality above. 
But, according to the completeness relation (|TT|), the difference involves the Lie derivative of 
i^£( A )/3, which is linear homogeneous in £(m and projects to zero.) 

Commutation relations for D A follow most easily from the definition fl2~4l) . For a scalar field 

/, 

2D [B D A] f = [Ci (B) ,C i(A) }f = C[t {B) 4 [A) ]f = e AB u a d a f, 

where we have recalled the well-known result that the commutator of two Lie derivatives is the 
Lie derivative of the commutator (i.e., Lie bracket), and made use of (Oil)- 
Consider next the operation on a 2- vector X a . We have from fl24|), 

D B D A X a = e { ; ) C e(B y {b) D A X b ) 

= 4% fl) {A^ (A) (e^)}. 

The projection tensor A^ can be replaced by 5@ , because the Lie derivative, operating on the 
difference, gives terms proportional to £^ or i^E)ai which project to zero, noting ([B2]). Thus, 

D B D A X a = e { p ] Ci (B) Ci iA) (e" a) X a ). 

We can now proceed exactly as for the scalar case to derive the commutator. The result (gen- 
eralized to an arbitrary 2-tensor) is 

[D B , D A }X a b \ = e AB C w dX\ ; . 
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In particular, 

[Dbi DA\g a b = 2€AB^(a;b) 

Recalling (B3), this may be written 

D[BK A ] ab = - €AB ^(a;b), 

which was used in (153). It contracts to 



1 

D[ B K A] = -e AB u . a . 

These last two identities also play a role in symmetrizing — or, more properly, recognizing 
the implicit symmetry of — the raw expressions for Rab and R Aa Bb that emerge from the Ricci 
commutation relations. The manifestly symmetric expressions listed in (p8|) and Sec. 13 have 
been symmetrized with the aid of these identities. 

To conclude, we note the rule for commuting Da and the two-dimensional covariant derivative 
V a . The commutator [D A , V a ], applied to any 2-tensor, is formed by a pattern similar to its 
two-dimensional covariant derivative, but with T bc replaced by 

d aKc = 2K A(b-,c) ~ K Abc ' a - 

As examples: 

[D A ,V a ]X b = X d D A T b da} 

[D A ,V a g bc ] = -2{D A T d a{b )g c)d = -2K Abc , a 

The justification for the rule is that the partial derivative d a (applied to any two-dimensional 
geometrical object) commutes with both d A and the two-dimensional Lie derivative L s d , sothat, 

by (0), 

[D A ,d a ] .V';;;. 0. 
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